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Abstract 

Using a total of 2.74 xlO''' decays of the 'il^{2S) collected with the CLEO-c detector, we present 
a study of Xcj -^ iV ■, where V = p^ lOJ-, (j). The transitions Xci -^ 7P^ and Xci -^ 1^ are observed 
with B{xc\ -^ 19^) = (2.43 ± 0.19 ± 0.22) x lO"^ and B{xc\ ^ 7^) = (8.3 ± 1.5 ± 1.2) x 10"^ 
In the Xc\ -^ IP^ transition, the final state meson is dominantly longitudinally polarized. Upper 
limits on the branching fractions of other XcJ states to light vector mesons are presented. 



Radiative decays of charmonium provide a rich context in which the interplay between 
theory and experiment can advance our understanding of quantum chromodynamics (QCD). 
The radiative decays of the J/ip that proceed through annihilation of the cc quarks are of 
particular interest for spectroscopy as they provide a gluon-rich hadronic system recoiling 
against the radiated photon. Such experimental channels are thought to be ideal for search- 
ing for bound states of gluons (glueballs); however, in order to interpret experimental data 
for these decays, one must have an understanding of radiative transitions of J /if) to light 
P-wave isoscalar (/j) states. In the case of the scalars (/o), the picture is complicated by 
the uncertainty in the structure and properties of the many observed experimental states. 
The radiative decays of P-wave charmonium {xcj) to light quark vector states (p°, cj, and 
(j)) provide an independent, complementary, cc-annihilation decay where the properties and 
structure of the final state hadronic system are well-known, which may be useful in validating 
theoretical techniques. 

In this Letter, we present the first observation of radiative decays of the xa to the light 
vector mesons p^ and uo. The measured rates for these decays are an order of magnitude 
higher than those predicted by Gao, Zhang, and Chao [Ij with perturbative QCD (pQCD) 
methods. 

The data used in this analysis were taken with the CLEO-c detector operating at the 
Cornell Electron Storage Ring (CESR) [2], which provided symmetric e^e~ collisions at the 
ip{2S) center-of-mass. The detector, described in detail elsewhere [21 I5J, features a solid 
angle coverage of 93% for charged and neutral particles. The charged particle tracking 
system operates in a 1.0 T axial magnetic field and achieves a momentum resolution of 
~0.6% at p = 1 GeV/c. The CsI(Tl) calorimeter attains photon energy resolutions of 2.2% 
aX E^ = 1 GeV and 5% at 100 MeV. Two particle identification systems, one based on 
ionization energy loss (dE/dx) in the drift chamber and the other a ring-imaging Cerenkov 
(RICH) detector, are used to identify pions, kaons, and protons. Detection efficiencies are 
determined using a GEANT-based |1] Monte Carlo (MC) detector simulation. 

To enhance photon energy resolution and reduce background, photon candidates are 
required to be detected in the barrel portion of the calorimeter (| cos6'| < 0.81) and must be 
spatially separated from the trajectories of charged tracks that have been extrapolated to the 
calorimeter. We form 7r° candidates from two photons whose invariant mass M (77) is less 
than three standard deviations from the nominal vr^ mass. For charged particles, we require 
a hit in at least 50% of the radial layers intercepted by the trajectory of the particle in the 
drift chamber, the x^/d.o.f. for the fit to the hits be less than 50, and the charged particle be 
consistent with originating from the e~^e~ interaction. To reduce backgrounds from Bhabha 
events, we additionally require that | cos 6*1 < 0.83 for reconstructed charged tracks. Defining 
ax as the number of standard deviations the measured dE/dx is away from the expected 
dE/dx for a particle of type X, we identify pions and kaons by requiring a,, < 4 and ax < 4, 
respectively. In addition, we utilize information from the RICH detector: Lx = — 21n£x5 
where Cx is the likelihood that the signature in the RICH is from a particle species X. For 
kaon candidates with p > 800 MeV/c that produce a signal in the RICH detector, we require 
Lk - L^ + crj^ - al <0 and Lk - Lp + a]^ - ap < 0. 

We reconstruct the exclusive decay 'ipi/lS) — > 7iXcj; Xcj — > IhV , where 71 (7h) designates 
the characteristic low (high) energy photon in the signal topology and V is either a p*^, a;, 
or (f) candidate. The p°, cj, and (f) candidates are reconstructed in the tt+tt^, Tr+Tr^Tr", and 
K^K~ decay modes, respectively. A four-constraint kinematic fit is performed to the entire 
event which forces the decay products to be consistent with the known four-momentum 



of the initial iIj{2S). Candidates that have x^/d-o.f. < 5 for this fit are retained. In the 
rare case that an event has more than one candidate, only the candidate with the smallest 
X^/d.o.f. is kept. The kinematically-fitted four-momenta of the decay products are used for 
subsequent analysis. To suppress multi-body hadronic decays of the ip{2S) we require that 
|M(7i7h) - M(7rO)| > ISMeV/c^ and |M(7i7h) - M(r/)| > 25MeV/c2. 

In the search for Xcj —>■ 7P° it is necessary to eliminate the copious background from 
e"'"e~ —>■ {e~^e~ or /i^/i~) where the leptons are misidentified as tt+tt" and radiated photons 
fake 7i and 7h. This background can be effectively eliminated by placing requirements 
on the opening angles of the two photon and pion candidates in the laboratory frame: 
—0.70 < cos6'^+,r- < 0.90 and |cos6'^j^,J < 0.98. To further suppress backgrounds from 
Bhabha events, the total detected energy in the calorimeter is required to be less than 
90% of the center of mass energy. An additional background arises from decays of the 
type ip{2S) -^ 7h^';^' -^ -yin^iT^ and is suppressed by requiring |M(7i7r+7r^) — M{ri')\ > 
15MeV/c2. 

Our general strategy for extracting the signal is to select events using the invariant 
mass of the candidate vector meson and then plot the transition photon energy E{'~f\) for 
events passing these selection criteria. The signal for Xco, Xa, and Xc2 decay will appear 
as peaks in -E(7i). The signal selection criteria for the three unique final states are 0.50 < 
M{7i+7i~) < 1.10GeV/c2 {xcj -^ 7P°), 0.75 < M{n+n-n'^) < 0.82GeV/c^ {xcj ^ 7^), and 
1.01 < M{K+K-) < 1.04GeV/c2 {xcj ^ 70)- 

The distribution of il>{2S) -^ 'JXcj transition photon energy -^(71) is shown in Figs. fl](a)- 
(c) for Xcj -^ iV, where V = p^,uj, and 0, respectively. Clear signals are observed for the 
Xci —>■ IP^ and Xci ~^ 1^ transitions. To extract the event yield from the spectra we first 
obtain a signal shape for each of the nine XcJ -^ lV transitions using an MC simulation 
of the signal where the mass and full width of the Xcj are taken from Ref. [5]. The MC 
simulation is subjected to the same kinematic fitting and analysis requirements as the data. 
Each of the distributions in Fig. [T] is fit to a linear background shape and a sum of three signal 
shapes, one for each of the Xcj states. The two parameters that describe the background and 
the normalization for each of the Xcj photon lines are allowed to float in the fit. The fitted 
yields are summarized in Table [Tj By examining the change in the fit likelihood when the 
signals yields are forced to zero, we estimate the significance of the Xci -^ 7P° and Xci -^ 1^ 
signals to be much greater than 5a, while the significance for Xci -^ 70 is less than 3a. 
These estimates do not include systematic uncertainties, discussed below, that may affect 
the significance of the yield. 

Our signal yield can be potentially biased by background from real Xcj decays, which peak 
in -^(71), that are partially reconstructed thereby faking our signal. Fortunately, hadronic 
decays of the type Xcj ~^ 7r° + (p°, u, or 0) are forbidden by C-parity conservation, otherwise 
they would certainly contribute a substantial peaking background to our signal. Other 
hadronic decays such as Xcj -^ K~^K~7t^ or Xcj -^ vr+vr" 71*^71° are allowed. In general, these 
either do not peak in vector meson invariant mass or require the loss of multiple neutral 
particles, and are consequently suppressed by the requirements placed on hadronic candidate 
invariant mass or x^/d.o.f. of the kinematic fit. In fact, using an MC simulation that models 
all ^{2S) and Xcj hadronic decays, we observe no such peaking backgrounds. Figure [2] 
shows the invariant mass distributions for p° and u: candidates in the xa region of -£^(71) 
- as is evident from the sideband regions, the bias due to non-p° or non-uj backgrounds 
is small. Nevertheless, for those channels where we have sufficient statistics to do so, we 
adopt a data-driven approach to estimate this bias. In the Xcj -^ 7P° case, we generate 
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FIG. 1: The ip{2S) -^ IXcJ transition photon (71) energy distribution for (a) XcJ -^ "fP^, (b) 
Xcj -^ 7^, and (c) Xcj — > 7'/' candidates. The data are shown by the points; the fit (described in 
the text) is shown as a sohd Une. The background component of the fit is indicated by the dashed 
Hne. 

background-subtracted -^(71) spectra by fitting the p° yield in bins of -^(71). Repeating this 
procedure with variations of the background parameterization in the p°-candidate invariant 
mass spectrum resulted in a maximum deviation from the nominal efficiency-corrected yield 
of —2% (—50%) for Xci(Xc2) -^ 7P°- The nominal analysis was also repeated while altering 
the selected region in M(7r^7r~). Changes in the efficiency-corrected yield for the Xci{Xc2) — *■ 
7P° signal ranged from -1% to +2% (-20% to +20%). For Xci -^ 1^, we extract the yield from 
a fit to the -£'(71) spectrum obtained by selecting events in the cu-candidate invariant mass 
sideband, 850 < M{n^n^n^) < 920 MeV/c^ (shown in Fig. ^, and conservatively assume 
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FIG. 2: Invariant mass of (a) p^ and (b) uj candidates for events that pass all analysis criteria 
and have £'(71) consistent with tpC^S) -^ JiXci (150 < £'(71) < 200 MeV). The points are data 
and the solid line is signal MC scaled to the yield extracted in the nominal fit. The signal region 
is indicated by the solid arrows. The sideband region for uj candidates (described in the text) is 
indicated by the dashed arrows. 



that this yield, 3.1 ± 3.2 events (8% of our signal yield), is equivalent to the background in 
the LV signal region in our nominal analysis. In addition, we repeat the analysis for various 
selected regions in M(7r+7r~7r°). In both cases, changes in the efficiency-corrected yields for 
Xci -^ 1^ were never larger than ±8%. In all cases described above, we find no statistically 
significant evidence for a bias in the efficiency-corrected yield. The central values for the 
(insignificant) biases are used as a quantitative estimate of our uncertainty, summarized in 
Table [Tj In all other channels we conservatively estimate the upper limit on the rates by 
assuming that all observed events are signal. 

The efficiency for each mode (see Table [l]) is obtained using an MC simulation that 
models the initial polarization of the ^(25') and the appropriate electric-dipole (El) angular 
distribution for the '0(25) -^ IXcJ transition photon. The decay Xcj —* iV is simulated 
uniformly in phase space except for the Xci -^ 7P° and Xci —^ 1^ decays. Here we modify 
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the MC to reflect the measured polarization, described in detail below. The efficiencies 
include branching fractions of the final state vector meson [5] and the detection efficiency 
for the initial transition photon 71. 

To obtain the product branching fractions B{ip{'^S) -^ IXcj) x ^iXcj -^ iV) we divide 
the yield by the product of the efficiency and the number of 'ip{2S) in our data sample, 
2.74 X 10^ [6]. The final Xcj —>■ jV branching fractions are obtained by dividing by the 
appropriate El transition rate, B{ip{2S) — > 7Xco) = 9.2%, B{ip{2S) -^ 7Xci) = 8.7%, or 
i3(-?/'(25') —^ 7Xc2) = 8.1% |;5J, and are summarized in Table [1} 

In addition to the impact of the biases described above, we explore several other sources 
of systematic uncertainty. Generous variations in the background parameterization used 
to fit the spectra in Fig. [I] produced variations no larger than 2%, 25%, and 2% for the 
Xci -^ 7P°5 Xc2 —>■ 7P°5 and Xci —>■ 1^ yields, respectively. Making significant changes in 
our event selection criteria produced variations in the efficiency corrected yield of 5% and 
8% for the xa — ^ 7P° and Xa — ^ 7^ channels, and we assign these respective systematic 
uncertainties to each of the Xcj -^ 7P° and Xcj — ^ 71-1^ rates. For Xcj — ^ 70, where we do 
not see a significant signal, we assume a systematic error due to event selection of 5%, the 
same as xa -^ IP^-, which has a similar topology. The errors in the track (photon) detection 
efficiency are assumed to be 1% (2%) per track (photon) and fully correlated across all 
tracks (photons). The number of ip{2S) in our data sample is known with 2% precision [6]. 
The uncertainty in polarization of the vector meson, assumed to be the maximum difference 
between phase space and either longitudinally or transversely polarized decays, introduces 
a 5%, 10%, and 8% error for the Xcj -^ 7(p°; ^, and 0) efficiencies, respectively, with the 
exception of Xci —^ 7P° and xa -^ 1^ modes where the polarization is measured and the 
resulting efficiency error due to uncertainty on this measurement is 1% and 3%, respectively. 
All of the ip{2S) —>■ JXcj rates have a relative 5% uncertainty |5j. The total systematic errors 
are summarized in Table [ij Upper limits are scaled by (1 + 5), where 6 is the total relative 
systematic error. 

The ratio of transverse (A = ±1) to longitudinal (A = 0) polarization of the vector meson 
A±/Aq can be measured by examining the distribution of events as a function of cos G, where 
is defined as the angle between the vector meson ffight direction in the Xcj rest frame 
and either the vr"'" direction in the p° rest frame or the normal to the decay plane in the 
u rest frame. Modulo detector acceptance, longitudinal (transverse) polarization exhibits a 
cos^ G (sin^ 0) dependence. The distributions of cos G are shown in Fig. pi where, for the 
Xci -^ 7P° case, the data (points) are obtained by ffiting the invariant mass spectrum of the 
vector meson candidate in bins of cos in order to eliminate potential contamination from 
non-p° decays. The Xci -^ 7^ candidates are plotted by requiring 150 < -E(7i) < 200 MeV. 
The individual transverse (dark gray) and longitudinal (light gray) components to which 
the data are fit are obtained from MC simulation, and the best ffi, floating A±/Ao and 
the overall normalization, is indicated by the total solid histogram. In principle, the decay 
amplitudes to the two polarization states can interfere; this interference is neglected in the 
fit. The fits give A±/Ao = Omst'^.'otl-o'oll for xa - 7P° and A±/Ao = OArt^.f+'oil 
for Xci —>■ 1^1 where the second, systematic error is obtained by assuming the estimated 
background contributes entirely to the longitudinal or transverse component. 

In summary, we present the first observation of radiative decays of Xci to light vector 
mesons. We find B{xci -^ 7P°) = (2.43 ± 0.19 ± 0.22) x 10"^ and B{xci ^ 7^^) = (8.3 ± 
1.5 ± 1.2) X 10~^. The measured rates are significantly higher than those predicted by a 
calculation using pQCD [IJ, for which the leading-order decay mechanism is annihilation 
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of the cc quarks into a light-quark pair that radiatively decays to ^yV. The longitudinally 
polarized structure of the Xci -^ 7P° decay parallels that measured in the decay of the 
corresponding light-quark axial- vector /i(1285) — *■ 7p° by VES [7J. This observation may 
suggest that the enhanced rate is due to the presence of a virtual light-quark axial-vector 
meson in the decay. The branching fraction measurements and upper limits presented in 
this Letter provide input to cross-check current and future calculations of radiative decays 
of charmonia that are important for spectroscopic interpretations of experimental data. 
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